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Abstract 



In the framework of general two-Higgs-doublet models, we calculate the branch- 
ing ratios of various inclusive charmless b decays by using the low energy effec- 
tive Hamiltonian including next-to-leading order QCD corrections, and examine 
the current status and the new physics effects on the determination of the charm 
QJ \ multiplicity n c and semileptonic branching ratio B$l- Within the considered pa- 

■ rameter space, the enhancement to the ratio BR(b — > sg) due to the charged-Higgs 

penguins can be as large as a factor of 8 (3) in the model III (II), while the ra- 
\ tio BR(b — ► no charm) can be increased from the standard model prediction of 

2.49% to 4.91% (2.99%) in the model III (II). Consequently, the value of B SL and 
n c can be decreased simultaneously in the model III. The central value of B$l will 
be lowered slightly by about 0.003, but the ratio n c can be reduced significantly 
from the theoretical prediction of n c = 1.28 ± 0.05 in the SM to n c = 1.23 ± 0.05, 
1.18 ± 0.05 for m H + = 200, 100 GeV, respectively. We find that the predicted n c 
and the measured n c now agree within roughly one standard deviation after taking 
into account the effects of gluonic charged Higgs penguins in the model III with a 
relatively light charged Higgs boson. 
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I. Introduction 



In the forthcoming years, experiments at SLAC and KEK B-factories, HERA-B and other 
high energy colliders will measure various branching ratios and CP- violating asymmetries 
of B decays Jl], ||. The expected large number of B decay events ( say 10 8 — 10 9 ) may 
allow us to explore the physics of CP violation, to determine the flavor parameters of 
the electroweak theory, and to probe for signals or evidences of new physics beyond the 
Standard Model (SM) [1-6]. 

Among various B meson decay modes, the decay b — > 37 and b — > sg have been, 
for example, the hot subject of many investigations J7|, since these decay modes may 
be affected by loop contributions from various new physics models. Great progress in 
both the theoretical calculation || and the experimental measurement enable one to 
constrain the new physics models, such as the two-Higgs-doublet model (2HDM) ||10|| , the 
minimal supersymmetric standard model [[LlJ and the Technicolor models |12j| . 

For many years, it appeared that the SM prediction for the semileptonic branching 
ratio B$l [0 is much larger than the values measured at T resonance and Z°-peak |I4 , 
[0| . More recently, the theoretical predictions have been refined by including full 0(a s ) 
QCD corrections [[16], [17|- These progress, consequently, have lowered the predicted Bsl 
and now adequately reproduce the experimental results ||15[ . However, the measurements 
of B,sl obtained at the T(4S) and Z° resonance are still disagree slightly ||18|| . Besides 
the Bsl problem, there is another so-called "missing charm puzzle" |15|, the charm 
multiplicity n c measured at CLEO and LEP [Tj| |2l|( especially at CLEO, the T resonance 
) is smaller than the theoretical prediction. Among various possible explanations for the 
missing charm/ Bsl problem, the most intriguing one would be an enhanced B — > X noc h ar m 
rate due to new physics beyond the SM Jl9| . An enhanced b — > sg can decrease the values 
of both n c and the Bsl simultaneously ||19| . The large branching ratio BR(B — > rj'X s ) 
reported recently by CLEO ]2TJ provided a new hint for enhanced b —>■ sg. Besides those 
explanations based on the SM p2| , new physics interpretation for this large ratio is also 
plausible [p3fl . 

In a previous paper |f24| , we calculated, from the first principle, the new contributions 
to inclusive charmless b quark decays b — > sg,b —>■ sqq from the gluonic charged-Higgs 
penguin diagrams in the so-called Model III: the two-Higgs-doublet model with flavor 
changing couplings P5L ^6] . In the considered parameter space, we found that the branch- 
ing ratio BR(b — > sg) (q 2 = 0) can be increased by roughly an order of magnitude, which 
is much larger than that in the ordinary 2HDM's P7|| . In plf , however, we used the 
language of form factors F\ and F2 and took into account the QCD corrections partially 
by using the a s (rrib) directly to calculate the branching ratios. 

In this paper, in the framework of general 2HDM's, we will calculate the branching 
ratios of various inclusive charmless b decays by using the low energy effective Hamilto- 
nian including next-to- leading order (NLO) QCD corrections ||, and investigate the new 
physics effects on the theoretical predictions for both Bsl and n c . 

This paper is organized as follows. In Sec. II, we describe the basic structures of the 
model III, extract out the Wilson coefficients, draw the constraint on parameter space of 
the model III from currently available data. In Sec. Ill, we calculate the branching ratios 
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BR(b — > sg) and BR(b — > g'gg) for q' £ d,s and g G w, cf, s in the model III and II with 
the inclusion of NLO QCD corrections. In Sec. IV, we examine the current status and new 
physics effects on the determination of B$l and n c . The conclusions and discussions are 
included in the final section. 



II. The general 2HDM's and experimental constraint 



The simplest extension of the SM is the so-called two-Higgs-doublet modelsflTUl. In such 



models, the tree level flavor changing neutral currents(FCNC's)are absent if one intro- 
duces an ad hoc discrete symmetry to constrain the 2HDM scalar potential and Yukawa 



Lagrangian. Lets consider a Yukawa Lagrangian of the form 1 26 



where (pi (i 
Qu. {Uj, R ) 



= 1, 2) are the two Higgs doublets of a two-Higgs-doublet model, 0i j2 = ir 2 (p\ 2 , 
with i = (1,2,3) are the left-handed isodoublet quarks (right-handed up- 



type quarks), Dj R are the right-handed isosinglet down- type quarks, while % f and £ 4 J° 
= 1,2,3 are family index ) are generally the nondiagonal matrices of the Yukawa 
coupling. By imposing the discrete symmetry 



-A, Ui 



(2) 



one obtains the so called Model I and Model II. In Model I the third and fourth term 
in eq. (|lj) will be dropped by the discrete symmetry, therefore, both the up- and down- 
type quarks get mass from Yukawa couplings to the same Higgs doublet (pi, while the (p 2 
has no Yukawa couplings to the quarks. For Model II, on the other hand, the first and 
fourth term in Eq.(fT]) will be dropped by imposing the discrete symmetry. Model II has, 
consequently the up- and down-type quarks getting mass from Yukawa couplings to two 
different scalar doublets (pi and (p2- 

During past years, the models I and II have been studied extensively in literature 
and tested experimentally, and the model II has been very popular since it is the building 
block of the minimal supersymmetric standard model. In this paper, we focus on the third 
type of 2HDM ||25|| , usually known as the model III p5| , f26fl . In the model III, no discrete 
symmetry is imposed and both up- and down-type quarks then may have diagonal and/or 
flavor changing couplings with (pi and (p2- As described in [2~B|, one can choose a suitable 
basis (H° , H 1 , H 2 , H^) to express two Higgs doublets J2(| 



^2\v + H° + ix 



^2 V H 1 + iH 2 



(3) 



and take their vacuum expectation values as the form 



< 



> 




< 02 >= 0, 



(4) 
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where v = (y/2Gp) l / 2 = 246GeV. The transformation relation between (H°,H 1 ,H 2 
and the mass eigenstates (H°,h°,A°) can be found in |p6 . The H ± are the physical 



charged Higgs boson, H° and h° are the physical CP-even neutral Higgs boson and the 
A is the physical CP-odd neutral Higgs boson. After the rotation of quark fields, the 



Yukawa Lagrangian of quarks are of the form |26 



Cy = 'hjQu.o.rj.R + TjfjQ i; L(piDj iR + C/Ij.oJ'j.h + ZijQi*<hDj, R + H.c, (5) 

where r^f correspond to the diagonal mass matrices of up- and down-type quarks, while 
the neutral and charged flavor changing couplings will be |26| [] 



C U,D _ \f m i m j x £U,D _ AJ,D PU _ AJ\r _ T/ tD , n 

Sij — ~ A ij} ^neutral ~ S > ^charged ~ S V CKM } ^charged ~ v CKMQ , \V) 



where Vckm is the Cabibbo-Kabayashi-Maskawa mixing matrix | 28| , i,j = (1,2,3) are 
the generation index. The coupling constants Ajj are free parameters to be determined 
by experiments, and they may also be complex. 

In the model II and assuming tan/3 = 1, the constraint on the mass of charged Higgs 
boson due to CLEO data of b -> 57 is M H+ > 350 ( 200 ) GeV at the LO (NLO) level 
29] , |3~0|1 . For the model I, however, the limit can be much weaker due to the possible 



destructive interference with the SM amplitude. 

For the model III, the situation is not as clear as the model II because there are more 
free parameters here. As pointed in PBJ, the data of K°—K° and B^—B^ mixing processes 
put severe constraint on the FC couplings involving the first generation of quarks. One 
therefore assume that, 

Kj = ^dj = 0, for j = 1,2,3 (7) 
Imposing the limit in Eq.([7|) and assuming all other Ay parameters are of order 1, 



Atwood et al. [31] found a very strong constraint of M H + > 600GeV by using the CLEO 



data of b — > 57 decay available in 1995. In Ref.[[$2[], Aliev et al. studied the b — > 57 decay 
in the model III by extending the NLO results of the model II |3(| to the case of model 
III, and found some constraints on the FC couplings. 



In a recent paper J33|, Chao et al., studied the decay b — > 57 by assuming that only 



the couplings \ u and Ay, are non-zero. They found that the constraint on M H + imposed 
by the CLEO data of b — > sj can be greatly relaxed by considering the phase effects of A« 
and Xbb- The constraints by B° — B° mixing, the neutron electric dipole moment (NEDM), 
the Z°-pole parameter p and Rf, give the following preferred scenario |33| : 



|A tt |<0.3, |A bb |^50, M A o w M h o = 80 - 120CeV; 80GeV < M H+ < 200GW (8) 

In the following sections, we will calculate the new physics contributions to the inclu- 
sive charmless decays of b quark in the Chao-Cheung-Keung (CCK) scenario of model III 



[33]|. Such model III has following advantages: 



1 We make the same ansatz on the ^ i -' D couplings as the Ref. p9|. For more details about the definition 



cU,D 



of £ u < n one can see Ref.|§. 
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1. Since we keep only the couplings X u and none zero, the neutral Higgs bosons do 
not contribute at tree level or one-loop level. The new contributions therefore come 
only from the charged Higgs penguin diagrams with the heavy internal top quark. 

2. The new operators Og 10 and all flipped chirality partners of operators Ox ... xo a $ 
defined in do not contribute to the decay b — > sy and other inclusive charmless 
decays under study in this paper. 

3. The free parameters in this model III are greatly reduced to Xtt, Xbb an d Mh+. 

In order to find more details about the correlations between M H + and couplings X t t,bb 
by imposing the new CLEO data of b — > sj, we recalculate the decay b — > sy in the model 
III. For the sake of simplicity, we do not consider the less interesting model I further in 
this paper. 

The effective Hamiltonian for B — > X s y at the scale /x = 0{m h ) is given by || 



U eff (b -> ay) = -^V t * s V tb 



Ci(ji)Qi(ji) + C 7T (/i)Q 77 + C 8G (/i)Q 8G 



(9) 



The explicit expressions of operators Qi_ 6 ,Q 77 and Qsg, as we U as the corresponding 
Wilson coefficients Ci(Mw) in the SM can be found for example in [|J. 

In the model III, the left-handed QED magnetic-penguin operator and the left- 
handed QCD magnetic-penguin operator Q$ G may also play an important role, 



Q77 = ^m^ll-^fe^, (10) 

Qsg = ^rn b s a a^{l- l5 )T^b p G% (11) 

In the SM and ordinary 2HDM's, both operators Q 7 and Qg G are absent because one 
usually assume that m s /mb ~ 0. In the model III, however, these two left-handed op- 
erators may contribute effectively because the Wilson coefficients and C$ G may be 
rather large to compensate for the suppression of m s /mb- 

In Ref. ^M, we calculated the b — > sg decay in the model III from the first principle 
and obtained the corresponding form factors F\ and F 2 . Following the standard procedure 
and using the Feynman rules in the model III |]26|| , we evaluate the Feynman diagrams 
for both b — > 57 and b — > sg decay as shown in Fig.|I|, extract out the Wilson coefficients 
Ci(M w ) at the energy scale M w by matching the full theory onto the effective theory, 

Q(M W ) = (i = 1,3,4,5,6), (12) 
C 2 (M W ) = 1, (13) 

771 

CU M w) = -^D(y t )\X tt \ 2 , (14) 
C*{M W ) = C 7l (M w ) SM - ^A(y t )\X a \ 2 + l -B(y t )\X tt X b b\e\ (15) 



12 w/l 1 2 



171 

c L m {Mw) = -j^- b D (yt)M> ( 16 ) 

Cf G {M w ) = C 8G {M W ) SM - ^D{y t )\X tt \ 2 + \E{y t )\X tt Xbb\e\ (17) 
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with 

C 71 {M W ) SM = (18) 

C 8G (M W ) SM = (19) 

where x t = mf/M^, y t = mf/Mj I+ , the phase angle 9 = 9b — 9 t , while 9b ( 9 t ) is the phase 
angle of Xbb (A#). When compared with the Eqs. (18,19) of Ref . |33| , the second and third 
terms in Eqs.(|l5|) and fll7D have an additional factor of 1/2, since £^- D used here has as 
additional factor l/y/2. The Inami-Lim functions (A, B, D, E) are of the form, 

7x - 5x 2 - 8x 3 2a; 2 - 3x 3 
12(1 -x) 3 + 2(1 -x) 

^ 4(1 -,)3 + 2TT^F lQsM ' (21) 

-3y + 5y 2 2y - 3y 2 

= wr^-6(r^F logM ' (22) 
= 4^-2a^F log ^ (23) 



^(*) = 10,1 ^ + o/i 34 10 SN, (20) 



The Wilson coefficients given in Eqs.(|T^-[l7|) contained the contributions from both the 
W^-penguin and if ^penguin diagrams. 

It is easy to see that both Cf 7 (il%) and C^ G (M W ) in Eqs.(|T|) and (H) will be doubly 
suppressed by the ratio m s /mb and \\tt\ 2 when |A tt | is small as preferred by the data 
of NEDM fl33f. For typical values of relevant parameters, say \X tt \ = 0.3, |A^| = 40, 



9 = 0° and M H+ = 200 GeV, One finds numerically that C^(M W ) « C% G (M W )) « 10" 5 , 
while C5,(M[y) ~ Cj^M^) ~ 0.8. Consequently, the left-handed Wilson coefficients are 
much smaller than their right-handed counterparts and therefore will be neglected in the 
following calculations. 

At the lower energy scale \i = 0(rrib), the Wilson coefficients Cj(/i) for the decay 
b — ► 57 at the leading order are of the form 

CM = E (J = !,•••, 6), (24) 

i=l 

C^Cai)^ = ^^(M^ + ^ryi-^C^M^ + EWS (25) 

^ i=l 

C 7 » m = V^C 7 R y (M w ) + l(r ] ^-r ] ^)C^ G (M w ) + j:h lV a % (26) 



i=l 



where 77 = a,(Jl%)/a,(//), and the scheme- independent numbers a,, fcjj and Ttj can be 
found in ||. 

Using the effective Hamiltonian, the branching ratio of b — > 37 at the leading order 
can be written as, 



B R(b -> S7 ) (m) = ^^-■^|C , 77 ( A i) m | 2 5i2(6- > oeP) J (27) 









K 6 


2 
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where \x = 0{m^), BR(b — > cez/) = (10.7 ± 0.4)% is the measured semileptonic branching 
ratio of b decay, and f(z) is the phase space factor, 

f(z) = l-8z 2 + 8z 6 - z 8 -2Az 4 \og[z] J (28) 

where z = m p c ole / 'm^ ole '. It is straightforward to write down the branching ratios BR(b — > 
57) for the SM and model II. 



In the numerical calculations, the following input parameters |15| |35| will be used 
implicitly: 

M w = SOAlGeV, M z = 91.187GeV, a em = 1/137, 
a s (M z ) = 0.118, G F = 1.16639 x 10- 5 (GeV)~ 2 , 
m s = O.UGeV, m c = lAGeV, m b = A.8GeV, 



(5)_ 
MS 

A = 0.84, A = 0.22, p = 0.20, rj = 0.34, (29) 



mt = m- t {m t ) = 168GeV, = 0.225, 



where A, A, p and rj are the Wolfenstein parameters of the CKM mixing matrix. m~i{mt) 
here refers to the running current top quark mass normalized at /i = m t and is obtained 
from the pole mass m v t ole = 176 GeV. For the running of a s , the two-loop formulae 
will be used. 

Fig.2 shows the branching ratios BR(b — > sj) in the SM and models II and III, 
assuming Aft = 0.3,Afcb = 35, 9 = 0°,30°, tan/3 = 1. The horizontal band between two 
dotted lines corresponds to the CLEO data |: 2x 10" 4 < BR(b -> 57) < 4.5 x 10~ 4 . 
The short-dashed line is the SM prediction, and the long-dashed and solid curve show the 
ratio in the model III for 9 = 0°, 30°, respectively. The dot-dashed curve shows the same 
ratio at the leading order in the model II. From the Fig.2, the lower and upper limit on 
M H + in the model III can be read out: 

\8hGeV < M H+ < 238GW, for 9 = 0°, 

215GeV < M H+ < 287GeV, for 9 = 30° (30) 

These limits are consistent with those given in Eq.(|8|). If we take into account the errors 
of theoretical predictions in model III, the corresponding mass limit will be relaxed by 
about 20 GeV. 

From above analysis, we get to know that for the model III the parameter space 
^ij — 0, for ij 7^ tt, or bb, 

|A tt | = 0.3, |A fe6 |=35, # = (0°-30°), M H+ = (200 ± 100)GeV, (31) 

are allowed by the available data. For the mass M H +, searches for pair production at 
LEP have excluded masses M H + < 77GeV Combining the direct and indirect limits 
together, we here conservatively consider a larger range of 100 GeV < Mjj+ < 300 GeV, 
while take = 200 GeV as the typical value. 
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III. Inclusive charmless b quark decays 



In this section, we will calculate the new physics contributions to the two-body and three- 
body inclusive charmless decays of b quark induced by the charged Higgs gluonic penguin 
diagrams in the models II and III. 

A. b — > s gluon decay 

The branching ratio of b — > sg at the leading order can be written as, 



BR{b -> sg) 



\v; s v tb 


2 &a e (p) 






2 


nf(z)K(z) 



with 



C 8G ^) SM 



C 8G (v)\ 2 BR(b^cev) 



rj^C 8G (M w ) SM + J2hiV ai 
i=i 

i=\ 



(32) 

(33) 
(34) 



where rj = a s (Mw) / a s (fi) with \i = 0(mb), and the numbers aj and hi can be found in [Q. 
The factor ^(z) contains the QCD correction to the semileptonic decay rate BR(b — > ceu) 
371 BR [39|]. To a good approximation the k(z) is given by K 



k(z) = 1 



2a a (/x) 
3tt 



7T 



+ 2. 



(35) 



And an exact analytic formula for k(z) can be found in ref . p8| . 

For b — > dg decay, one simply substitutes V t * s by V t * d in Eq. (p2[) . For the model II, one 
simply replaces Csg(^) in Eq.(|32]) with Cg G as given in [pTf . 

Fig. 3 shows the branching ratios of BR(b — > sg) in the SM and the models II and 
III, assuming \ tt = 0.3, A&& = 35, and 6 = 0°, 30°. The dots line in Fig. 3 is the SM 
prediction BR(b — > sg) = 0.27%, while the short-dashed curve shows the branching ratio 
BR(b sg) = 0.81% in the model II assuming tan/5 = 2 and = 200 GeV. In 

the model III, the enhancement to the ratio BR(b — > sg) can be as large as an order 
of magnitude: £#(& -> sg) w 2.34%, 4.84% for M H + = 200, 100 GeV respectively, 
as illustrated by the long-dashed and solid curves in Fig.3. The model III is clearly 
more promising than the model II to provide a large enhancement to the decay b — > sg. 
Although the current enhancement is still smaller than ~ 10% as expected, for example 
in Refs.|l9|, j2"3"| , such a significant increase is obviously very helpful for us to provide a 
reasonable solution for the problems such as the " missing charm puzzle" or the deficit 
Bsl, as being discussed below. 



B. Three-body charmless b quark decays 

Within the SM, the three-body inclusive charmless b quark decays have been calculated 
at LO and NLO level for example in refs.|^|, |23|, ^(J. In Ref.||, Lenz et al. took into 
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account the NLO QCD corrections from the gluonic penguin diagrams with insertions of 
Q 2 and the diagrams involving the interference of the Q$g with Qi-6 0- 

The standard theoretical frame to calculate the decays b — *> sqq for q G {u, d, s} is 
based on the effective HamiltonianJ]]], 



H eff (\AB\ = 1) = ^ \Y, C j (v c Q c 3 + v u Q]) - v t 



CjQj + C 8 Qsg 

j=3 



+ h.c, (36) 



where v q = V*V q b and the corresponding operator basis reads: 



with q = u and q = c, and 



Qi = (s a qp) V - A (qpb a ) V - A , (37) 
Q2 = {s a q a )v-A{qpbp) V -A, (38) 



(39) 
(40) 
(41) 
(42) 
(43) 



Q 3 = 


(s a b a 


q=u,d,s 


Q 4 = 


{Sab/3, 


V-A X! (QpQa)v-A, 
q=u,d,s 


Q 5 = 


(s a b a 


W-A {Wp)v+A, 
q=u,d,s 


Q 6 = 


(s a bp] 


V-A Y (Wcdv+A, 
q=u,d,s 


Q&G — 


9s 
8tt 2 


m b s a a^(l + l5 )T^bpG 



where the Qi and Q2 are current-current operators, the Q3 — Qq are QCD penguin oper- 
ators, while the Qsg is the chromo-magnetic dipole operator. 

For the SM part, we will use the formulae presented in || directly. For the new physics 
part in the models II and III under study here, we take into account the new contributions 
from charged-Higgs gluonic penguins by using the Wilson coefficient Csg^) 111 as given 
in Eq.(|34"D in the calculation, this coefficient comprises both the SM and the new physics 
contributions. All other Wilson coefficients remain unmodified. 

When the NLO QCD corrections are included, one usually expand the decay width to 
order a s , 

T(b^sqq) = T^ + ^-(AT cc + AT peng + AT w + AT 8 )+0(a 2 s ), (44) 

where denotes the decay rate at the LO level, while the second part represents the 
NLO QCD corrections. We here use the renormalization-scheme(RS) independent terms 
Ar cc , Ar pen9 and Ar^. For the convenience of the reader, the explicit expressions of 
Ar cc , Ar pe „ 9 and Ar^y will be given in Appendix. The term Ar§ in Eq.([l4]) ( which will 
be defined below in Eq. (Fj^) ) is already RS independent^, |33|. For the three-body decays 
b — > dqq one simply substitutes s by d in Eqs. (|36|-|4"4D . 
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At the NLO, the RS dependent Wilson coefficients Cj(fi) are given byfl35 



6. 



(45) 



where Cf^ are the RS independent LO Wilson coefficients, and Cf^ are the RS dependent 



NLO corrections [35] 



cfM = E^ ai 



i=3 



cfV&) = J2l e jiVE {xt) + f j i + g j iV]v at 



(46) 
(47) 



i=3 



where r\ = a s (M w )/a s (fib), x t = mf/Myy, the function E Q (x t ) and all the numbers a iy 
kjh e jii fjii an d Qji can be found in |p5| . The NLO QCD correction Cj is RS dependent 
and can be split into two parts: 



•,6. 



(4* 



k=l 



where parameters Jjk are usually RS dependent, Cjipp 1 ' is RS independent, and the 
precise definitions of the terms in Eq.(^) can be found for example in [pi] . The terms 
involving Jj k will be absorbed into Ar cc and AT peng to make the latter scheme indepen- 
dent. 

In the leading order the decays b — > sss, sdd, dss and ddd are penguin- induced pro- 
cesses proceeding via Qss and Qsg, while 6 — > dim and 6 — > sww also receive contributions 
from Qi and (^2- Combining both cases, the decay width at the LO level can be written 

as m 



r (0) 



i,i=3 



-2t E Q OJ q 0J ReK<)6ii 



i=l,2 
J=3,...6 



(49) 



with £ = 1 for q = u and t = for q — d,s. The coefficients read 



1 fi7T 3 

= ^ r _ / c /$ 3 (2vr) 4 (g j )(0)(g j )(0)' 



m 



6 



.7' 



(50) 



with Qi t 2 = Q12 here. Setting the final state quark masses to zero one finds | 

b; 



J 1 + r/3 for i,j < 4, and i + j even , 
1 1/3 + r for i,j < 4, and i + j odd , 

bm = 1 , ^56 = &65 = 1/3. 



(51) 
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Here r = 1 for the decays b — > ddd and b — > sss, in which the final state contains two 
identical particles, and r = otherwise. The remaining fc^'s are zero. 

Now we turn to study the contributions from the interference of the tree diagram with 
Q 8 with operators Qis, as shown in Fig.3 of Ref.@. The tree-level correction Ar 8 is 
already at the order of a s and is given by 



AT* 



G 2 F m\ 
32vr 3 



Re 



-tvlv t C 8G (n) 



in 



Y.cf% + \v t \ 2 c, G ^) III Y.cf ) b j s 

i=3 



in the model III, where Csg(^) has been given in Eq.([34|) with \x 



(52) 

0(m&). For the 



case of the SM and model II, simply replace CgG((ib) in with the appropriate Csg (/■*&)■ 
The definitions and numerical values of coefficients bj$ can be found in ||. As men- 
tioned previously, the Wilson coefficient C^q now comprises the contributions from both 
the W-penguin and the charged- Higgs penguin diagrams. In this way, the new physics 
contributions are taken into account. 

For the b quark decay rates one usually normalize them to the semileptonic decay rate 
of b quark, 

qWl) r(6 



ql 



r(6 
r(6- 



cev e ) 
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Q 9) 
F(b — > ceu e ) ' 

r(6 -> sgg) 



919293 



T(b^ceu e ) ' ' syy r(6->ceP e )' (53) 

for the sake of eliminating the factor of mf common to all b decay rates. One also define 
the charmless decay rate of b quark as 

(54) 



r 4 



^ ] (j"dgg T ' sqq) T ' sg ~\~ 1 dg T sgg 27" ue ~\~ 7\ 
q=u,d,s 



where rare radiative decays, for example b 
semileptonic decay rate takes the form 



r(b 



cev P 



sj, have been neglected. To order ac s , the 
'"' b] V ch \ 2 f{z)K{z) (55) 



G 2 F m 5 

1927T 3 



where the factors f(z) and k(z) have been given in Eqs. (|28f) and Q35I) . 

To calculate we also need explicit expressions of r. 
finds [0, 

'\ - z 2v ' 
1 + k{z) - k(0) + 6 1 



nr.) ^sgi f "dg and fsgg- For T ue 



one 



\v ub 


2 1 


\v cb 


2 f{z) 



/(*) 



1 



A, 



mi 



(56) 



where A2 = 0.12GeV 2 encodes the chromomagnetic interaction of the b quark with light 
degrees of freedom, and the factors of f(z) and k{z) have been given in Eqs. fl28D and (j35]). 
From Eq.(B3), we get 



r dg 



\v t :v tb \ 2 




\v cb \ 2 


7rf(z)K(z) 


\v td v tb \ 2 


8a s (//) 


\Vcb\ 2 


7t/(z)k(z) 



\c 8G (ri\ 2 , 

\CSG^)\ 2 - 



(57) 
(58) 
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For r sgg} we use the formulae as given in [40. 24 

1 3« s (/i) 2 



r sgg 



|K & | 2 lGir 2 f(z)K(z) 



£ VWMxW)?, (59) 



i=u,c,t 



where X{ = m 2 /M^, the functions fi(xi,q) can be found for example in |2"4| . In the 
numerical calculation, we assume that q 2 = ml/2. Since the new contribution to the 



decay b — ► sgg due to the charged Higgs penguin is negligibly small [2H], we do not 
consider the new physics corrections to this decay here. In Ref.|]], the authors did not 
include r sgg in the estimation of ry. We here will include this mode, since its branching 
ratio is rather large |40 , 24 1 , as shown in the Table [I]. 



The corresponding branching ratios for two-body and three-body charmless b decays 
are defined as 

BR(b->X) = r x ■ BR(b — > ceu e ) exp } (60) 

where ratios rx have been defined previously. In the numerical calculations, BR{b — > 
ci/ e ) exp = 10.70% will be used Q. 

By using the input parameters as given in Eq.fl2"§D and assuming \ \ tt \ = 0.3, |Aw>| = 35, 
Mh+ = 200GeV and 9 = 0° or 30°, we find the numerical results of the decay rates and 
the branching ratios for various charmless b quark decays and collect them in Table |l|. 
We also show the corresponding results in the model II assuming Mh+ = 200 GeV and 
tan/? = 2. For larger tan/3 the new physics contributions in model II will become smaller. 
ABR in Table [1| is defined as 



ABR(b -+X)= [BR(b -> X) - BR(b -> X) bM \ /BR(b -> SM) (61) 

Fig.^ shows the mass dependence of the branching ratios BR(b — > s qq) with g e 
{u, d, s} in the SM and model III, using the input parameters in Eq.(p9|) and assuming 
\\ tt \ = 0.3, | | = 35, and 9 = 30°. In Fig.f|, the three curves ( horizontal lines) are 
the theoretical predictions in the model III ( SM ) for q = u, d, s, respectively. For 
Mh+ = 200GeV, as listed in Table [I], the enhancement to the decay mode b —>■ duu 
is only 4.7%, but the enhancements to other five three-body b quark decay modes are 
rather large: from ~ 30% to ~ 70%. In the model II, however, the new contributions 
are negative and will decrease the branching ratios slightly, from —0.3% to —13.5% for 
different decay modes. 

Fig.|5] shows the branching ratio BR(b — > no charm) in the SM and models II and III, 
using the input parameters in Eq.(^Sj) and assuming \X tt \ = 0.3, | A&5 1 = 35, and 9 = 0°, 
30°. The dots line in Fig.|5] is the SM prediction BR(b — > no charm) = 2.49%. The 
short-dashed curve shows the the ratio in the model II, BR{b — > no charm) = 2.98% 
(3.23%) for M H + = 200 (100) GeV and tan/? = 2. The long-dashed and solid curve show 
the theoretical predictions in the model III: BR(b — > no charm) = 4.67% (4.91%) for 

2 For more details, one can see the discussions about the semileptonic branching ratios of b decay in 
next section. 
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Table 1: The rates r and branching ratios in the SM and models II and III, assuming 
|A«| = 0.3, |A 66 | = 35, M H+ = 200GeV, tan/? = 2, and 9 = 0° or 30°(the numbers in 
parenthesis). We also use BR(B — > X c ev e ) exp = 10.70%. as given in Eq. (|63"D . 



SM Model III Model II 



decay mode 


r 


BK(7o) 


r 




ABK{7o) 


ABK(7o) 


b - 


-> duu 


0.051 


0.545 


0.052 


0.554 


1.6 


-0.3 










(0.053) 


(0.571) 


(4.7) 




b - 


_ 

-> ddd 


0.0005 


0.00b 


0.00078 


0.103 


68.2 


1 O A 

-13.4 










(0.0007) 


(0.010) 


(59.0) 




b- 


->■ dss 


U.UUUo 


U.UUo 


U.UUUyo 


U.UUo 


fie i 
oo. ( 


-lo.o 










(0.0009) 


(0.008) 


(59.5) 




b- 


-> suu 


0.018 


0.192 


0.027 


0.286 


49.0 


-9.6 










(0.0237) 


(0.255) 


(32.8) 




b- 


^ sdd 


0.019 


0.206 


0.030 


0.322 


56.0 


-11.0 










(0.0285) 


(0.307) 


(48.5) 




b- 


^ s ss 


0.016 


0.168 


0.024 


0.262 


56.7 


-9.9 










(0.0232) 


(0.250) 


(49.2) 




b- 


-> s# 


0.025 


0.270 


0.192 


2.065 


663.6 


202.3 










(0.217) 


(2.339) 


(765.0) 




b- 


-> dp 


0.00092 


0.010 


0.007 


0.070 


663.6 


202.3 










(0.008) 


(0.086) 


(765.0) 




b- 


-»• s ££ 


0.070 


0.757 


0.070 


0.757 






b- 


-> it ez/ e 


0.013 


0.144 


0.013 


0.144 






b- 




0.013 


0.144 


0.013 


0.144 






b- 


-> it rz/ T 


0.004 


0.0004 


0.004 


0.0004 






b- 


-> no charm 


0.23 


2.49 


0.43 


4.67 


87.6 


20.3 










(0.46) 


(4.91) 


(97.3) 
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M H+ = 200 GeV and = 0°, 30°, respectively. For the model III with M H+ = 100 GeV, 
one finds that BR(b — > no charm) = 7.27% (7.60%) for 9 = 0°, 30°, respectively. 

It is easy to see from Fig.[5| and Table |1] that the new physics enhancement to the 
branching ratios of three-body charmless b quark decays in the model III is much larger 
than that in model II within the parameter space considered. 



IV. n c and B$l 

The ratio Bsl is the average over weakly-decaying hadrons containing one b quark. For 
the CLEO experiments running on the T(4S) resonance, the average is over B + and 
B and their charge conjugate hadrons. For the experiments running on Z° resonance, 
however, the average is over B + , B°, B® and N^. 

The charm multiplicity n c is the average over the b-hadrons produced in the given 
environment. CLEO and LEP collaborations presented new measurements of inclusive 
b — ► c transitions that can be used to extract n c . One naively expect n c = 1.15 with 
the additional 15% coming from the tree-level decay chain b — > uW~ — > ucs. This 
expectation can be verified experimentally by adding all inclusive b — > c branching ratios, 
and counting twice for the decay modes with 2 charm quarks in the final state. 

In this section, we will investigate the new physics contributions, induced by the 
charged Higgs penguins in the models II and III, to the ratio Bsl and the charm multi- 
plicity n c . 



A. n c and Bsl- experimental measurements 

The Bsl deficit was first point out in around 1994 J05] when the theoretical prediction 
was considered to be difficult to produce B S l < 12% while the 1995 CLEO data on T(4S) 
resonance was B S l = (10.49 ±0.46)% [pp . In the following, we use the 1998 Particle Data 



Group value [15 



B SL = (10.45 ±0.21)% (62) 



as the measured Bsl on T(4S'). 



For the experiments on the Z -peak, all the four LEP collaborations [43, 44, 45, 46 



reported their measured values of the ratio Bsl as listed in Table |[ The seventh row 
shows the averaged result of the ratio B S l on Z°-peak0: B b SL = (10.66±0.17)%. This B b SL 
on the Z°-peak can be converted to T(4S*) value by multiplying a factor of tb/t^ = 1.026: 
Bsl = (10.94 ± 0.19)% (Z° corrected). In fact, there is still a 2a discrepancy in ratio 
Bsl between the high energy Z° value and the low energy T(45) value. The average of 
the Z° and T(4S*) values of B S l is 

B SL = (10.70 ± 0.21)% Overall average (63) 

where we conservatively chose 0.21 as the overall error of the measured Bsl- 
Nb is in turn the mixture of Ab(udb), Eh(ws6), 'E.b(dsb) and flb(ssb). 

4 We here made an arithmetic average over four results as done in but the newest L3 data Q 
has been used here in the average. 
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Table 2: Recent CLEO and LEP measurements of the ratio Bsl- 



B S l (%) 






Experiment 


10.45 ±0.21 






T(45) PDG98 |lj 


11.01 ± 0.10(stat.) 
10.65 ± 0.07(stat.) 
10.16 ±0.13(stat.) 
10.83 ±0.10(stat.) 


±0.30{syst] 
±0.25{syst.] 
±0.30(syst] 
±Q.2olsyst] 


+° f 2 (model) 
±° f 3 (model) 


ALEPH 95 p| 
DELPHI 99 |44] 
L3 99 |45| 
OPAL 99 g§ 


10.66 ±0.17 






Z°-peak 


10.94 ±0.19 






Z° corrected 


10.70 ±0.21 






overall average 



As for the charm counting, the value of n c measured at the T(4S f ) f2(J is still smaller 
than that measured at Z°-peak [15] : 



_ f 1.10 ±0.05, T(4S), 
c " \ 1.20 ±0.07, Z°-peak l ° J 

The average of the T(4S) and Z° result leads to 

n c = 1.14 ± 0.04, (Z° + T(4S)) (65) 

B. n c and Bsl> theoretical predictions 

Within the SM, the basis of the prediction for B$l and n c is the assumption of quark- 
hadron duality. The estimation for various inclusive decay rates is usually performed by 
using the heavy-quark expansion(HQE) |47j and the perturbative QCD in the framework 
of operator product expansion. The HQE allows to relate the inclusive decay rate of B 
meson to that of the underlying b quark decay process: T(B — > X) = T(b — > x)+0(l/mf). 

The theoretical prediction for Bsl with the inclusion of the 0(a s ) QCD corrections 
and the hadronic corrections to the free quark decay of order l/m| is currently available 



16, 171. The Bsl and n c can be defined as El n 



Bsl = — — , (66) 

J2i r d + r cUd + r c5s + r 4 

n c = 1 + - r J^Zll (67) 

Y.i r d + r cU d + r cSs + r 4 

where r ce = r CAt = 1, r CT = 0.25, and r cU d (r c5s ) is the rate of the decay mode b — > cud' 
(b — > ccs') where d' (s') is the appropriate Cabibbo mixture of d and s quarks. 
The has been defined and calculated in last section. In the SM, we have 

r 4 = 0.23 ± 0.08, (68) 
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where the error mainly comes from the uncertainties of the scale fi and the mass ratio 
m c /m b @. 

As is well known, the main difficulty in calculating B$l and n c is in the non-leptonic 
branching ratios r C ud and r c5s - For r C ud, a complete NLO calculation has been performed 
[ |T6] which gives 

r cUd = 4.0 ± 0.4, (69) 

where the error mainly comes from the uncertainties of the scale /z, the quark masses m c 
and m b and the assumption of quark-hadron duality [16]. Furthermore, the error of the 



estimation for r c£s is generally considered to be larger than that for r cU d- The enhancement 



of b — > ccs due to large QCD corrections is about 30% |I6| . Such enhancement will 
decrease the value of B SL , but increase the size of n c . 

Using the on-mass-shell sheme, the SM theoretical predictions for B$l and n c at the 
NLO level are 

B SL = (12.0 ±0.7 ±0.5 ±0.2+^)%, (70) 
n c = 1.24 =f 0.05 ±0.01, (71) 



as given in Ref.|T6|0; and 



_ / (12.0 ±1.0)% fi = m b , 
* SL ~ \ (10.9 ±1.0)% n = m h /2, U j 

/ 1.20 =F 0.06 fi = m b , 
Hc \ 1.21 =f 0.06 n = m b /2, 1 } 



as given in Ref.||17|| with the error mainly result from the variation of the scale /i and 
m c /m b . 

Comparing the observed and predicted values of B SL and n c , one can see that: (a) after 
considering all the corrections, the theoretical values of B$l now come down and more or 
less consistent with the measurement, but unfortunately at the expense of boosting n c ; (b) 
the central value of n c in Ref . [jl6[| is higher than that in Ref.[l7j], although two predictions 
are agree within errors; (c) there is still 2.8 a discrepancy between the n c measured by 
CLEO and the theoretical prediction 0: 1.10 ± 0.05 against 1.24 ± 0.05 



If we'd like to drop down the large uncertainty in the calculation for b — > ccs' decay 
mode, we can eliminate the ratio r c5s i from the expression of B$l and n c and find, 



2 - (2.25 + r cUd + 2r 4 ) B SL (74) 



which is a linear correlation between B SL and n c . Using the values for B SL (p3|) , r cUd (p9|), 
and (]68|), one finds 

n c = 1.28 ±0.05, (75) 



5 The last and largest error of Bsl comes from the uncertainty of the renormalization scale \x\ while 
the main error of n c is the the uncertainty in mi, |lq|. 
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for Bsl — (10.70 ± 0.21)%. The overall uncertainty of this prediction of n c should be 
smaller than that as given in Eqs.([H]) and ([731). The 2.6 a discrepancy between the n c in 



Eq.flTSP and n c measured at T(4S") motivated proposals of new physics which will enhance 
and in turn decrease n c . That is what we try to do here. 

As shown in Table [I], the ratio will be increased significantly after taking the new 
physics effects into account, which will in turn decrease both Bsl and n c accordingly. 
From Eq.flTID and using the values for B SL flB3|), r cUd (|5T3|), and Ta (0), one finds 

f 1.23 ± 0.05 for M H + = 200GeV . , 

nc ~\ 1.18 ±0.05 for M H + = lOOGeV [ ' 

for 9 = 30° and \x = mj. The fi- and ^-dependence of n c is rather weak: the central 
value of n c will go down (up) by only ~ 0.01 for /i = m b /2 {6 = 0°). For B SL in the 
model III, the agreement between the prediction and the data will be improved slightly 
by a decrease 0.003 (0.005) for Mh+ = 200 (100) GeV due to the inclusion of new physics 
contributions. In the model II, the resulted decrease for n c (Bsl) is only 0.01 (0.001) 
and plays no real role. Most importantly, one can see from Eqs. (|6^j65| , [75|j76|) that the 
predicted n c and the measured n c now agree within roughly one standard deviation after 
taking into account the effects of gluonic charged Higgs penguins in the model III with a 
relatively light charged Higgs boson, as illustrated in Fig|| 



V. Summary and discussions 

In the framework of the general two-Higgs doublet models, we calculated the charged- 
Higgs penguin contributions to (a) the rare radiative decay b — > S7; (b) the inclusive 
charmless decays b — > q'g and b —* q' qq with q' 6 {d, s} and q G {u, d, s}; (c) the charm 
multiplicity n c and semileptonic branching ratio Bsl- 

In section II, we studied the experimental constraint on the model III from the CLEO 



data of b — » 57 decay. With the help of previous works [|26|, [U], [32], [33|, we found the 
parameter space of the model III allowed by the available data, as shown in Eq. (PT|) . 

In section III, we firstly calculated the new physics contributions to the decay b —>■ sg 
and found that the branching ratio BR(b sg) can be greatly enhanced from the SM 
prediction of 0.27% to 2.34% (4.84%) in the model III for M H + = 200 (100) GeV, as 
illustrated in Fig. 3. Such a significant enhancement is clearly very helpful to resolve the 
missing charm/ Bsl problem appeared in B experiments. 

Following the method of Ref . , we then calculated the new physics contributions to 
three-body inclusive charmless decays of b quark due to the interference between the oper- 
ators Qi-q and Q$g- The Wilson coefficient C%q in Eq. (|34D now describe the contributions 
from both the W ± and QCD penguins, the latter is the new physics part we focus in 
here. From numerical calculations, we found that: (a) the new physics enhancement to 
the decay b — > duu is only ~ 1.6% since this mode is dominated by the tree diagrams; (b) 
the branching ratios of other five three-body b decay modes are strongly enhanced by the 
new charged Higgs penguins: 30% to 70% increase can be achieved within the considered 
parameter space. The new contributions to the corresponding branching ratios in the 
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model II is, however, small in size and negative in sign against the theoretical predictions 
in the SM. As shown in Table [l] and Fig.|5|, the ratio BR(b — > no charm) can be increased 
from the SM prediction BR(b — *> no charm) = 2.49% to BR(b — > no charm) = 4.91% 
(7.60%) in the model III for M H+ = 200 (100) GeV. 

In section IV, we studied the current status about the theoretical predictions and 
experimental measurements for the semileptonic branching ratio of B meson decay B$l 
and the charm multiplicity n c , and calculated the new physics contributions, induced by 
the charged Higgs penguins in the model III (II), to both Bsl and n c . With an enhanced 
ratio BR{b — > no charm), both the B SL and n c will be decreased accordingly: (a) the 
central value of B SL can be decreased slightly by 0.003 (0.005) for M H + = 200 (100) GeV; 
(b) the value of n c can be lowered significantly from the prediction n c = 1.28 ± 0.05 in 
the SM to n c = 1.23 ± 0.05, 1.18 ± 0.05 for M H+ = 200, 100 GeV, respectively. 

In short, the predicted n c and the measured n c now agree within roughly one standard 
deviation after taking into account the effects of gluonic charged Higgs penguins in the 
model III with a relatively light charged Higgs boson, while the agreement between the 
theoretical prediction and the data for Bsl can also be improved by inclusion of new 
physics effects. 
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Appendix: RS independent Ar cc , Ar peng and AlV 

For the convenience of the reader, we here present the explicit expressions of the RS 
independent NLO corrections Ar cc , AT peng and Ar^. For more details one can see the 
original paper 0. 

The term Ar cc in Eq.(f44"D describes the current- current type corrections proportional 



to c$.cQ0 



Glm 5 



2 



k=l 



(77) 



with t — 1 for q = u and t — for q — d,s, and the coefficients h%j and Jk% can be found 
in§. 

The term AT peng in Eq . (fill) describes the effect of penguin diagrams involving Q\ 2 

I. 



- G 2 F m\ 
peng = 1^" 



t £ Cfcf v u [v* cgij (x c )+v* ugij (0)} 
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i=l,2 
i=3,...6 



+ Re 



E cfcf j w 6 ifc + 1&| 2 E d^dfijub 



<jh 



fe=3,...6 



i=l,2 
j,fc=3,...6 



(7f 



with t = 1 for q = u and £ = for q = d, s. The explicit expressions of coefficients gij and 
Jki can be found in ||. 

Finally, ATw is given by 



Ar 



32tt 3 



i E W 



^(0) t^(i) 

* 3 



% + E H 1 

i,i=3 



6 ij 



* E 

i=l,2 
J=3,...,6 



rr(o) i r»(o) 



ite (» t ) 6. 



(79) 



where t — 1 for q = u and t = for q — d,s, the 6^ have been given in Eq. (Ell) . 
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Figure 1: The Feynman diagrams for the decays b — > s^y and 6 - 
2HDM's. The internal quarks are the upper type u, c and t quarks. 



sg in the SM and 




Q * 

100 150 200 250 300 

M H (GeV) 



Figure 2: Plots of the branching ratio BR(b — > 37) versus M#+ in the SM and models 
II and III. The short-dashed line is the SM prediction, and the band between two 
dots lines refers to the CLEO data. The dot-dashed curve shows the ratio in the 
model II, while the long-dashed and solid curve show the ratios in the model III for 
9 = 0°, 30°, respectively. 
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Figure 3: Plots of the branching ratio BR(b — > sg) versus M#+ in the SM and models 
II and III. The dots line is the SM prediction, the short-dashed curve shows the the 
ratio in the model II, and the long-dashed and solid curve show the ratios in the 
model III for 9 = 0°, 30°, respectively. 




Figure 4: Plots of branching ratio BR(b — > sqq) versus M H + in model III. The three 
curves ( horizontal lines) are the theoretical predictions in the model III (SM) for q = 
u,d,s, respectively. 
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Figure 5: Plots of the branching ratios BR(b — > no charm) versus M H + in the SM 
and models II and III. The dots line is the SM prediction, the short-dashed curve 
shows the the ratio in the model II, and the long-dashed and solid curve show the 
theoretical predictions in the model III for 6 = 0°, 30°, respectively. 
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Figure 6: Plots of Charm multiplicity n c versus M H + in the SM and model III for 
Bsl = 10.70%. The short-dashed line is the SM prediction, and the band refers to 
the data of n c = 1.14 ±0.04. The solid curve, the upper and lower dot-dashed curves 
together show the central value and the 1 a error of the theoretical prediction for n c 
in the model III. 
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